Effective tight-binding model for the iron vacancy ordered Kj^Fei.6Se2 
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We investigate the electronic structure of the ternary iron selenide Ki,Fei.6Se2 by considering the 
spatial symmetry of the y/E x \/5 vacancy ordered structure. Based on three orbitals of t2g, which 
are believed to play major physics in iron-based superconductors, an effective two-dimensional tight 
binding Ifamiltonian is constructed with the vacancy ordered structure being explicitly included. It 
is shown that the constructed band model, when combined with generalized Hubbard interactions, 
yields a spin susceptibility which exhibits both the block-checkerboard antiferromagnetism instabil- 
ity and the stripe antiferromagnetism instability. In particular, for large Hund's rule couplings, the 
block-checkerboard antiferromagnetism wins over the stripe antiferromagnetism, in agreement with 
the observation in experiments. We argue that such a model with correct symmetry and Fermi sur- 
face structures should be the starting point to model Kj;Fei.6Se2. The spin fluctuations at q=(7r,7r) 
suggest that interblock fluctuations of spins might play an important role in the mechanism of 
superconductivity occurring in this system. 

PACS numbers: 74.70.Xa, 74.20.Pq,74.20.Mn 



I. INTRODUCTION 

After the discovery of superconductivity in LaFeAsO-^ , 
the finding of iron-chalcogenide superconductor a-FeSe^;^ 
has stimulated another intensive studies on the iron- 
based superconductors. Although the iron-chalcogenides 
and iron-pnictides have similar crystal and electronic 

[ structures, they show many differences either in super- 
conductivity or in magnetism. For instance, the former 
has lower (~ 8K)'^ and a larger magnetic moment of 
the Fe ion in FeTei_a:Sea; '^fJ-B)^ while the latter, e.g. 
La-1111, has much higher Tc = 26K— but has a smaller 
magnetic moment 0.36/1^) for the Fe ion^. Even the 
magnetic orders are different, one is bi-collinear and the 
other is collinear. Recently, K, Cs or Rb intercalated 
FeSe superconductors Aj^Fe2-a;Se2 are found. It is shown 
that Tc of this system can be enhanced above 30K-^. In 
this system, as the atomic ratio of Fe:Se is not 1:1, which 
used to be in iron-chalcogenides, iron deficiency is pro- 
duced. As a result, in addition to the enhancement of 
Tc, a y/b X \/5 iron vacancy ordered pattern illustrated 

I in Fig. [T] is formed when the composition is close to 
KyFei.eSeai"— . The vacancy ordering is followed by the 
magnetic transition at lower temperature Tjv — 560K 
to the block antiferromagnetic (AFM) phase with a big 
moment S.Slyit^. Although the AFM state is observed 
in superconducting K0.8Fe1.6Se2, evidence of nanoscale 
phase separation from X-ray diffraction is reported'^. 

Previous theoretical works^^^^'^ on superconductiv- 
ity in the AyFei.6Se2 system were based on the band 
structure of KFe2Se2 in which the hole pocket around 
r is absent and only electron pockets are present at 
M. Such Fermi surface tomography was supported by 
ARPESi^""— . However, since there is no experimental 
evidencei^ii^ showing the existence of the -v/S x a/5 pat- 
tern in KFe2Se2, the validity of these approaches is ques- 
tionable. In fact, because of the -s/S x -\/5 pattern, the 



symmetry group changes from H/mmm to H/m and 
both the unit cell and the Brillouin zone (BZ) change 
as well. The underlying band structure should be very 
different. Indeed, the first-principles calculations of the 
-\/5 X \/5 vacancy ordered lattice structure^Si^i indicate 
that a hole pocket at F appears in the nonmagnetic state 
and as expected it loses the reflection symmetry in the x-y 
plane. The presence of a hole pocket at F point indicates 
that the physics that drives superconductivity could be 
very different. It thus calls for a close examination based 
on an appropriate Hamiltonian to model the ternary iron 
selenide KyFei.6Se2. 



In this paper, we construct a two-dimensional tight- 
binging model with three t2g orbitals for the K0.8Fe1.6Se2 
system with the -y/S x \/5 vacancy ordered pattern be- 
ing explicitly included. Based on the general tight- 
binding model Ht with symmetry imposed by the va- 
cancy order, we fit the dispersion relation of Ht to that 
of the non-magnetic state from the first-principles cal- 
culations of Ref^S. Two hole pockets and two electron 
pockets emerges in the fitted tight-binding model. The 
constructed band model, when combined with general- 
ized Hubbard interactions, yields a spin susceptibility 
that shows a large peak around the q- vector (tt, tt) for 
the block-checkerboard AFM state. Furthermore, com- 
petition of different magnetic states is found but the 
block-checkerboard AFM state gets enhanced with larger 
Hund's rule coupling and wins out at the end. The im- 
plication of our results to the mechanism of supercon- 
ductivity is discussed. In particular, we argue that the 
spin fluctuations at q=(7r, tt) suggest that in analogy to 
spin-fluctuations in high Tc cuprates, the interblock fluc- 
tuations of spins might play an important role in the 
mechanism of superconductivity occurring in this system. 
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FIG. 1: (Color online) The schematic representation of the 
vacancy ordered lattice structure of the iron plane. 
The black dots denote Fe atoms and the dashed lines enclose 
the unit cells, x and y are the primitive vectors. 



II. THEORETICAL MODEL 

We start by considering the vacancy-ordered structure 
with right-handed chirality2£ as shown in Fig. [TJ In 
this ordered state, the unit cell includes four iron atoms 
(without considering Se), which we denote as I = A, B, 
C, and D and three t2g orbitals {d^z, dy^, and d-^y) are 
considered in each iron. Therefore we have 12 species 
of electrons in a unit cell. We will suppress spin indices 
and denote the electron operators collectively as a vec- 
tor hy 1p = (■0i,V'2: V's) ■with Ipr = idrA,drB,drC,drD), 

where drM denotes the annihilation operator of electron 
for orbital r at site M with r = 1, 2, 3 standing for d^z, 
dyz, dxy, respectively. We will use x and y as the coor- 
dinates of this system and x and y as the nearest Fe-Fe 
directions. 

In previous works, the DFT band structure of the par- 
ent compound KFe2Se2 in which only electron pockets 
appear at Mii"— is employed. However, the parent com- 
pound KFe2Se2 does not have the same hopping parame- 
ters as those in the vacancy-ordered Kj^Fei.6Se2. For ex- 
ample, there is no hopping between Fe atoms and vacan- 
cies. Here we shall strictly enforce the -\/5 x -\/5 vacancy- 
ordered structure and construct a tight-binding model 
with nearest neighbor (NN) and next-nearest neighbor 
(NNN) Fe-Fe hoppings, which will be classified as intra- 
and inter-cell ones. As the vacancy ordering appears, the 
refiection symmetry is lost but the four-fold-rotational 
symmetry is left intact. The system is invariant under 
90° rotations around the center of a unit cell, which is 
the position of Se (if the vacancy position is taken as 
the rotation center, the rotation has to be followed by a 
Pz operation: the reflection z — )■ —z). A general hop- 
ping Hamiltonian Ht with vacancy order being included 
can be written down by imposing the 90° right-handed 
rotation symmetry: di — > d2^ ^2 ^ ^'^i, ^3 ^ —d^ ac- 
companied with A^B,B^C,C^D,D^A. Due 
to its massive form, the general form of Ht is relayed 
to the Appendix. After Fourier transformation, in the 
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FIG. 2: (Color online) The band structure of the three-orbtial 
model with the \fZ x -\/5 lattice structure in the un-folded BZ 
(a) and in the folded BZ (b), and its correspoinding Fermi 
surfaces in (c) and (d), respectively. We have shifted the dis- 
persion (/i=0.557 in our model), such that £'=0 corresponds 
to the Fermi level. The bands shown are only those near the 
Fermi level, from the eighth to the eleventh (from low energy 
to high energy). 



momentum fc space, it takes the form 

i/i = ^^t(k)M(k)v.(k), 



(1) 



where V(k) = (V-iCk), '02(k), V-sW) with -(/^^(k) = 
(rfrA.k, f^rS.k, t^rCk, c?r_D,k) and M(k) is a 12 X 12 ma- 
trix. Detailed characterization of all hopping parameters 
are tabulated in TABLE [TT] and TABLE HH] in Appendix. 
These parameters are obtained by fitting energy disper- 
sions (in the folded BZ) to the results of X. W. Yan et 
al?'^ obtained by the generalized gradient approximation 
(GGA) in which the main features are four hole pockets 
at r and four electron pockets at X in the nonmagnetic 
state. Our fitting gives two hole and two electron pock- 
ets, which capture basic features of this system. Fig. [5] 
shows our fitting results in the unfolded BZ ((a) and (c)), 
and in the folded BZ ((b) and (d)). In the folded BZ, 
there are two hole pockets around (0, 0) and two electron 
pockets around (tt, 0); in the unfolded coordinate, one 
hole pocket will move to (7r,7r) and electron pockets to 
±(7r/2,±7r/2). 

As for the particle number, the stoichiometric com- 
pound Ao.8^'ei.6S'e2 gives Fe^"*". In other words, there 
are six electrons for each iron. Previous three-band model 
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TLB — no 


1/4 E/^/ 


m 


1.92 


1.74 


1.83 


n2 


1.74 


1.92 


1.83 


n-j, 


0.84 


0.84 


0.84 




4.50 


4.50 


4.50 



TABLE I: Particle number per Fe for different orbitals (r — 
1, 2, 3) and sites (7 = A, B, C, D). Due to the four-fold rota- 
tion symmetry, some numbers are equal. 



alytic continuity, the susceptibility becomes 
1 



ab.cd/ 



E A,^(k)A:^(k)^,(k + q)A3,(k + q) 



OJ - 



(4) 



We now include the effect of electron-electron interaction 
by considering the generalized Hubbard model, in which 
all interactions are on the same Fe atom. 



works for iron-pnictides claimed four electrons per Fe in 
the undoped state. 

In the GGA calculation^, the number of electrons en- 
closed by Fermi surfaces is about 0.642 electrons/cell, 
while the number of holes enclosed by Fermi surfaces is 
about 0.529 holes/cell. In our model at the symmetry 
point, the electron number per iron is 4.5 and hence the 
total number of electrons is 18 per cell. The particle den- 
sity for each orbital and site is listed in Table HI As we 
expect, due to symmetry, ni at site ^ or C (i? or D) is 
the same as n2 at site B or D {A or C), while is uni- 
form at every site. In addition, we found that there are 
about 0.52 electrons/cell and 0.52 holes/cell enclosed by 
Fermi surfaces. These numbers are close to those found 
in the GGA calculation. We note in passing that it is pos- 
sible to change the chemical potential and hopping scales 
so that the model is away from the symmetry point and 
numbers of electrons/holes per cell are closer to those 
obtained by the GGA calculation. However, since we do 
not find significant changes of magnetic properties, we 
shall be focusing on the symmetry point. 



III. MAGNETIC AND CHARGE RESPONSES 

Using the tight-binding model with the fitted param- 
eters found in the last section, we can analyze linear 
responses of the system. We shall first calculate the 
generalized susceptibility in the absence of the electron- 
electron interaction defined by 



ab.cd/ -y-v 



dre 



(^i,(q-^)^cd(-q-o))o (2) 



= J2 5ca(k, iw„)g6d(k-l-q, ia;„-|-iri„) 

Here the generalized spin operators are defined by S'^j, = 

V'l.t^'b.i and S;:^ = (S+^y = V'rf4'0c,t with the subscript 
(a, b, c, d) being the 12 orbital indices for electrons, and 
the Green's function gat is given by 



^ iujn - Ef^{k) 



(3) 



where is the band index and Aaf^ is the orbital-band 
transformation matrix, V'oC*) = 'I2^aii0<-)j f_i(^) ■ By an- 



i I=A,B.,C.,D [ a=l,2,3 

Jh 



a,b(a>b) 



U' ^ ) nal.inbl,i — 2JH^aI,i ■ ^bl. 



+Jc {dlJ^^^dlJ^^■^dbI,^idbI,i^ + h.c}j } 



(5) 



Here we simply use the same set of parameters for every 
site and orbital. Within this model, we calculate the 
random-phase approximation (RPA) susceptibilities for 
spin and charge 



1 - ^sXoi^,^) 

xo(q,^) 
i + rcXo(q.'^)' 



(6) 
(7) 



The vertices for spin sector are r^'^''^'^ — U, T^^ •^'^ — 
U', r7'^'^' = Jh, F^''^'^ = Jc, and for charge sector 

Yrr,TT ^ prr'.rr' ^ ^ 2Jh, F^'^'^' = 2U' - 

Jh, F^"^ ^ — Jc, where nonvanishing vertices are only 
between the same Fe, and r denotes orbitals and t ^ t' . 
In the following, we shall take the relations U' ~ U — 2Jh 
and Jc = Jh- 

In Fig. |3l we show the total DC susceptibilities per 
cell (four iron) defined by x'°*'''(q,0) = Es,tX''''"(q, 0). 
Here U—1.2eV and Jh = 0-2U are used. The black solid 
line is for the bare susceptibility Xo°*"'' the blue dashed 
line is the spin susceptibility, and the red dotted line is for 
the charge susceptibility. As expected, electron-electron 
interaction strongly enhances spin susceptibility ^'o*"' 
and induces a peak around (tt, tt). The Stoner instability 
for Jh — 0-2U is found to happen at ?7=1.5eV and such 
divergence of x*°*"' at (7r,7r) will result in the checker- 
board AFM pattern as experiments observed. Therefore, 
the fitted tight-binding Hamiltonian explains the experi- 
mental observations. On the other hand, the charge sus- 
ceptibility 1^*°*"' is not important here and is smaller than 
the bare one, which is consistent with results of Refi^. 

Next we investigate the effect of the Hund's-Rule cou- 
pling. As shown in Fig. SJa), for small U (take J7=0.5eV 
as a nominal example), values of show a mono- 

tonic behavior. However, in the large U case as shown 
in Fig. |li;b) ([/=!. 2eV), values of x*°*°' exhibit non- 
monotonic behavior. In particular, the shoulder around 
(tt, 0) at large Jh becomes a hump at Jh = 0. The hump 
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FIG. 3: (Color online) Total DC susceptibilities, black solid 
line: bare, blue dashed line: spin, red dotted line: charge. 
The interaction parameters are U = 1.2el/, Jh = 0.2?7, and 
U' = U "^Jh. 



at (tt, 0) indicates that there is a magnetic instability for 
striped AFM at about ?7=1.3eV, in competition with the 
block-checkerboard AFM. 

To further check the magnetic instability, we employ 
the Stoner criterion. In the multi-orbital system, the 
susceptibility is a matrix and magnetic instability is de- 
termined by the corresponding eigenvalues. The Stoner 
criterion requires one to find the first eigenvalue, A^, that 
reaches zero, i.e., As(q) = 0+, where A5(q) is the minimal 
eigenvalue of the inverse of X,,i^p^(q, 0) {Xs]ipa^'^^ 0) = 
X^^(q, w) — Ts] and q will be the magnetic ordering vec- 
tor. Fig. Ul^c) andm^d) show the behavior of As(q) ver- 
sus U . It is seen that at Jh — 0, shown in Fig. SJc), 
the first eigenvalue that touches zero occurs at q =(7r, 0). 
Hence the stripe AFM is the resulting magnetic phase at 
= 0, in consistent with our previous conclusion. At 
larger J//, as shown in Fig. Hfd) {Jh = 0.2[/), the mag- 
netic instability occurs at q =(7r,7r). From Fig. and 
mjd), we also find that the critical value of U , Uc, when 
magnetic instability occurs, depends on Jh as well. For 
Jh = 0, we find that Uc ~1.3eV, while for Jh = 0.2[/, 
we get Uc ~1.5eV. These results all suggest that large 
Hund's rule coupling stabilizes the checkerboard AFM 
state. 

We note in passing that in the above, we do not try to 
distinguish whether the magnetic instability occurs ex- 
actly at (tt, tt) or not. All of these magnetic states are 
classified as the checkerboard AFM state. In fact, be- 
cause the parameters adopted in Fig. |3] are for the sys- 
tem at the symmetry point, the magnetic instability does 
not happen exactly at (7r,7r). By changing the chemical 
potential, the wave vector of the magnetic instability can 
be shifted to be exactly at (tt, tt). This implies that the 
exact wave vector for the magnetic instability will gener- 



ally depend on the doping level of the system. 



IV. SUMMARY AND DISCUSSION 

In summary, in contrast to perturbative treatment 
of vacancies^^, we have constructed an effective tight- 
binding model for the Kj,Fei.6Se2 system by including 
exact symmetries of the Fe vacancy ordering structure. 
The tight-binding model includes three orbitals (d^z, dyz, 
and dxy), which are considered to be the most important 
orbits in iron-pnictides and iron-chalcogenides. Although 
this system shows a large moment^ and could be better 
described by including some localized moments, a proper 
tight-binding band structure is still required since iron- 
based superconductors so far are regarded as an interme- 
diate coupling system instead of being a strong coupling 
system. For example, recent experimental findings from 
thermal transport of Ka;Fe2-j^Se2 indicated it a weakly or 
intermediately correlated system=4. From these aspects, 
it is clear that our model captured the essential low en- 
ergy physics: two hole pockets around F and two electron 
pockets around X and Y in the folded BZ. Furthermore, 
the constructed band model, when combined with gen- 
eralized Hubbard interactions, yields a spin susceptibil- 
ity which exhibits both the block-checkerboard antiferro- 
magnetism instability and the stripe antiferromagnetism 
instability. In particular, for large Hund's rule couplings, 
the block-checkerboard antiferromagnetism wins over the 
stripe antiferromagnetism, in agreement with recent ob- 
servations in experiments. 

While so far in this work we only consider the magnetic 
instability of the ternary iron selenide system, our find- 
ings also provide some insight into possible mechanism 
for superconductivity occurring in this system. In partic- 
ular, the strong spin fluctuations at q=(7r, tt) could result 
in inter hole-pocket and inter electron-pocket (in opposite 
momenta) scatterings, which may lead to pairing with 
totally different symmetries of pairing. In real space, it 
implies that inter-block fluctuations of spins might play 
a similar role in analogous to spin-fluctuations in high- 
Tc cuprates. While our model has not yet accounted 
for superconductivity observed in this system, the fitted 
tight-binding model shall serve as a useful starting point 
for developing the correct theory. 
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FIG. 4: (Color online) Upper two panels: total spin susceptibility Xs?j{pa(Qi 0) different Jh in 17=0. 5eV (a), and in U—1.2eV 
(b). Lower two: Ai,(q), the minimal eigenvalue of the inverse of xI°r'pa{^j 0) three q vectors in Jh = (c), and in Jh = 0.2U 
(d). The magnetic transition happens when Aa(q) — 0+. 



Appendix A: The effective tight-binding 
Hamiltonian 

In this appendix, we will include details for construc- 
tion of the tight binding Hamiltonian. Following the 
symmetry argument given in the context and neglect the 
tetramer lattice distortion^", the tight-binding Hamilto- 
nian with NN and NNN hoppings can be written as 

Ht = H, + Hr2 + H3 + Hi2,3- (Al) 

Here is the on-site energy. H12 characterizes hopping 
among orbitals: d^z and dyz, while H3 is the hopping 



term for d^y and -ffi2,3 describes the hopping between 
d-xz/ dyz and dxy- 

We shall suppress the spin index for simplicity. To con- 
sider the effect of Se atoms above and below the Fe plane 
periodically, the transformation dsi^i — > (— l)'*'rf3/,i is 
included implicitly to make the hopping integrals site- 
independent. Due to symmetries imposed by the a/S x 
vacancy ordered structure, we find that the on-site ener- 
gies for d^ and dy^ are different and their difference will 
be denoted by A, while the on-site energy of d^ will be 
denoted by e. The on-site energy can be written as 



I 

i 

-(1 ^ 2)] (A2) 



6 



To describe hopping terms, we will adopt the notation 
t -B for intra-cell hoppings and t' for inter-cell hop- 

mn,R t-t- to m.n,R 

pings. The subscript mn are the orbital indices and R is 
the Fe-Fe direction. We note that because of the absence 



of reflection symmetry, the NN Fe-Fe hopping between 
dxz and dy^ is allowable now. By including all possi- 
ble terms allowed by symmetries, hopping terms can be 
generally expressed as 



Hl2 = 



i 



\2 — 

i 

idlA,i + C^lB,jC^lC,i + dlj^^d2B,i + C^2C,j'^2£),i^ 
+ tll,x (d\g j^diDA-x + dl^^d2C,i+y^ + t'li-y {d\^,^diCA+y + d\g^^d2DA-i 

tll,x+y {d\j^^dic,i + d\g^^d2D,i^ + tll,x-y {^XsAdlL 
+ ^'ll,x+y {d\A,id\. 



+ A\,x+y {d^\A,id\D,i+y + d^-^Q ^^d\B ,i-y + 4B,i<^2A,i-x + 4r>,i'^2C,i-|-x) 
+ A\,x-y {d^\A,id\B,i+x + C^Icit^lO.i-x + dls,id2C,i+y + 4 

+ ti2,x (dl^ id2B,i +d\^ ' - .+ , N 



^2B,id2C,i+v + dl 

. - ^ . Ac,id'i.D,i - dljj idici - dljj 

+ tl2,y (d\^^diBA + 4' 



4yl.j^^2C,i) 
^2D,id2A,i-y^ 



+ 



idiA,i^ 



-..,y yzA^i^^^^.^ ' Aci'^^Da - d\g^^d2C.i - d\jjj^d2A,?j 
' il2,x {dlB,id2D,i-x + dljj ^d2B,z+x - dly^ idiCA+y - 4, 
- tl2,x+y (d\A,id2C,z + d\c id2A,i - d^^^idlD^i - d^jj ^i 
' tl2,x+y (d\A,id2D,i+y + d\(. -d2B,i-y - d\g ^dlA,^-x 



d2C,i + d\cj_d2A,i - d\B,idiD,i - d^jj ^d 

(d^ 

,x+y 

idlD,i 



;,i -r 'J.xc,i"--iA,i - "2B,i"l£',i - "'2D,idlB,i^ 
d2D,i+y + d\(-,.d2B,i-y — d\g^dlA,i-x — d\jjidlC,i+x^ 
+il2,x-y {dlA idiD,i+y + d\(. ^d\B,i-y - d\g ^d2A,i-x - d\jj id2C,i+^ + /i-c.| , 



(A3) 



i 

i,x {dlA idiCi+y + dlB ^dsD,i-a ^ 

'' 't ' , .+ . ■-d\(j ,^diB,i-y + d\jj ^dic,i+^ +h.c)^ , 



+*33, 



^ , . „ , ^) + *33,x-Fy {dlA,idiC,i + d\i^ ^dsD,i^ 

,x+y {dlA^id3D,i+y + dlB id3A,i-x + die/" " ' '^^ - • ■ " 



(A4) 



andi/i2,3 = ^ |^13,^ (^3A,z^l^.i ~ ^SB.i^'^C.i — dl^ .diD,i + (^SD^i^'^^i) 

i 

+ ti3,y (dl^ ,^diD,t - d\g^^d2A,i - dlcidiB,! + 4n,i'^2C,i) 
+ i23,x {d\A,id3D,i — d\Q ^d3A,i — d\c^id3B,i + d\JJ ^d^c,i^ 
+ t2i,y {d\A idiB,i - dls idaci - rflc.i'^SD.i + 4D,iC^3A,i) 

+ *'l3,x (^D.i^SB.i+x - ^A.i'^SC.i+j/ - d\^ ^d3D,i-x + 4c,i'^3A,i-3/) 

+ ^13,1/ {d\^ ,^d3C,i+y - d\g .-d^D^i-x - d\c^id3A,i-y + 4z5,j'^3B,j+a;) 

+ il3,x+j/ (4c,i'^3A,i - d\j^ ^d3B,i - 4A,i^3C,j + 4B,i'^3D,i) 

+ *13,x-j/ (4B,i^3D,j - d\(-, ^dzA,i — d\jj ,id3B,i + 4A,j^3C,i) 

+ *'l3,x+y (^A.i'^SD.i+j/ - 4B,i^3A,i-x - 4c,i^3S,j-j/ + 4n,i'^3C,i+x) 

+ il3,x-y (4yl,j^3B,i+x - a'lc,,'^3i3,4-2; - d\g^^d3C,i+y + 4_D,i'^3A.i-2/) 

+ ^23,3;+i/ {d\A.idiDA+y ~ d\^^^d3B,i~y + d\g ^d-iA,i-x " d\^ ^d3C.i+x^ 

+*23,x-j/ (4A,i<^3B,i+x - d\Q ^dzD,i-x + <^lB,i<^3C,i+y - C^lD,i<^3A,i-j/) + /l-cj . 

I 



After Fourier transformation, the Hamiltonian is writ- 
ten in a matrix form as 



M(k) is given by 



i/t =^V^(k)M(k)^(k) 



(A6) 



M(k) 



where the basis vector is defined as before, 
V'(k) = (V'i(k),V2(k),V3(k)) with = 
(rfrA.k) c?TB,k) c?TC,k, rfrD.k) and the 12 X 12 matrix 



Mn(k) Mi2(k) Mi3(k) 
M|2(k) M22(k) M23(k) 

Mt3(k) M^yk) M33(k) 
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(A5) 



(A7) 



with elements being give by 



Mn(k) 



A 

2 



iii,x + i'li^^-yC * —-J + ^'ii,x+y^* 

(-11, y -|- i-u^x+y*^ 2 



til 



.x+y + ^'ll,a^ 



tll,x-y + t'u^x^ 



tll,x + A\,x-y^ 



ikx 



tll,y + t'll^x+y^ 



hi. 



x-y 



iii.x + t'n,x-ye^ 

A 



(A8) 



M22(k) 



_ A 
2 

^11,17 + tii,x+i7'^ 
' + *'ll,x^ 



ill, 

iii.x + t\ 



ll,x-y'^ 



-ik, 
-ik 



tll,y + til_3; 



A 
2 



'■^ tii,x + i'n.^-^e 

tll,x+y + ^ll,y^ 



ikx 



tii,x-y + i'ii,se*''« iii,x + i'l 1,^-176''^'' 
iii,x+|7 + i'n.i/e"'''" 



^ll.x + tii,x-17 

~"2 *11.17 + *'ll,x+17^ 

(-11,2/ i- ''11,5+ye 2 



(A9) 



M33(k) 



h3,x + tss.x+yC^ 

t33,x+y + ti^s^e-'f^y 



*33,x + ^S.x+ye''^^ i33,x+y + i33,xe " *33,x + ^S.x+ye'*^'' 



^33,3; + 43.x 



i33,x + ^SS.x+l/e' " t33,x+i? + 43,x« 



-ikx 



*33,x + 43,x+ye ""-^ 



(AlO) 



. i33,x + i33,x+y^ ' " t33,x+y + i'33,x^'' " *33,x + i33,x+y6 
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Mi2(k) = 



''12.x ''12,x+y^ 



-tl2,x + t[2,x+y'^ 



. -h2,y + tl2,x-y<^ -^12,^+1? + t[2^x'^ 



tl2,x+y - t'i2,x^ 
-tl2,y + A2,x-y' 



'-12,1/ ''12,x-y'^ 



tl2,x + tl2,x+y^ 



-tl2.x+y + t'i2.x 
tl2.x — tl2,x+y 





■e 

—ik 



(All) 



Mi3(k) = 





tl3,x + t-y^i xA 



'13, X -r 1-23,5+17 
13,3^+17 
tl3,y 



tl3,x+y ^13,y^ 



f' e 

''23,x-y^ 



t23,y + t'lZ.x-y^''^" 


h3,x - t'l3,x+y^^'^ " 

—ik-u f T T L fD^kx 

''13,x-y ^ ''13,x^ 



tl3,x+y + i'i3 " t23.x + t[ 



'13,1/ 



"'-13,2; ''23,x+17'^ 



^23,17 



13,Ic+i7"- 
~ ^13,S^ 
^13,11-1/^ 



ihj; 



-ikx 



(A12) 



and 



Af23(k) 



— ^23,^ — ^13,^+17^ 
~'tl3,x-y + ^13,^6 
^23,17 + t'lS.x-y':^ 



ik^ 



'.ky 



tl3,y + t23,x-y 



~^13,K ^ ^23,x+y 



-tl3,x+y + ^13,176'''" ^23,3; + t'l3,x+y^''^" 



''13,x-y '-13,5^ '-13,2; ^ ''23,5+17*^ 

'•13,2;+i; ''13,17'^ 

tl3,y 



-t 



23,17 



4-f f^'f'ky + _ _ 

' ''13,1e-i7'^ '•13,2;+i; 



' '•23,5-17^ 





(A13) 





!r = ^ 


R = y 


R^x+yR 


^x-y 


mn = 11 


-0.14 


-0.09 


0.03 


0.03 


mn = 33 


-0.05 




0.3 




mn = 12 













mn — 13 


-0.25 





-0.05 


0.15 


mn = 23 





-0.1 







Our fitting values are A=0.2 and e=0.55, and those of 
the hopping integrals are hsted in TABLE|lT]and TABLE 

ml 



TABLE IL Fitted intra-cell hopping parameters between NN 
and NNN. 



mn,H 


R = x 


R=yR^x+yR 


^x-y 


mn — 11 


-0.028 


-0.04 0.024 


0.024 


mn = 33 


-0.05 


0.35 




mn = 12 





-0.03 


-0.06 


mn = 13 


-0.375 


-0.05 


0.1 


mn = 23 




0.15 


0.05 



TABLE IIL Fitted inter-cell hopping parameters between NN 
and NNN. 
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